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Abstract 

We study how attractive boron correlations in boron-doped diamond affect the 
superconducting critical temperature. The critical temperature is obtained from 
the McMillan formula for strong coupling superconductors with the density of 
states evaluated in the dynamical cluster approximation. Numerical results for 
the cluster of 2 x 2 x 2 atoms show that attractive correlations lower the density 
of states at the Fermi level. We argue that this might explain experimentally 
observed differences in critical temperatures of 100 and 111 oriented films. 
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1. Introduction 

Diamond is a typical insulator with a band gap of ~ 5.5 eV. It has been 
intensively studied for a long time due to its compatibility with human tissues 
and unique physical properties like high transparency for a visible light or high 
thermal conductivity which are promising for future applications. When doped 
with boron, a p-type conductivity appears. As the doping concentration exceeds 
~ 4.5 x 10 21 cm -3 the system becomes metallic and eventually superconducting, 
with an unexpectedly high critical temperature on the order of a few Kelvin. 
[H, [H, HI The critical temperature depends on the method by which the sample 
was grown. 

Samples are most often prepared in thin layers, using the Microwave Plasma- 
assisted Chemical Vapour Deposition (MPCVD) method with growth orienta- 
tions 100 and 111; see [1] and [f| respectively. The 111 samples have higher T c 
than 100 samples. There are also bulk samples prepared with the High-Pressure 
High- Temperature (HPHT) method [f| with T c comparable to the 100 samples. 
Some experimental data are collected in Fig. [5] below. 

While the concentration of impurities is well under control in any of the 
above preparation methods, a concentration of binary correlations between 
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boron atoms is not and may cause the difference between samples. Here we 
discuss an effect of these correlations on T c . 

Our study is motivated by previous experimental studies. From the Raman 
scattering interpreted with the help of first-principle calculations Bourgeois et 
al suggested that boron atoms in boron-doped diamond form correlated multi- 
boron complexes, primarily dimers. Q This is confirmed by Mukuda et al who 
observed a strong NMR line of isolated boron impurities in 111 samples, while 
this line is rather weak in 100 samples. Q 

We will assume that boron atoms tend to form isolated dimers but not 
large clusters. This assumption is supported by ab initio investigations which 
show that the nearest-neighbour dimer is the most favourable of all two-boron 
configurations while an added third boron tends to remain unassociated with 
the dimer [l,[loj]. 

A dimer can be viewed as a molecule with bonding and anti-bonding states 
made of symmetric and anti-symmetric combinations of bound states of com- 
posing boron atoms. The symmetric state is deeper in the band gap than the 
boron state, theqrefore it does not contribute to N , the density of states at the 
Fermi energy, [ill 0,13 The energy of anti-symmetric states is above the isolated 
boron level and also does not contribute to Nq. The formation of dimers thus 
reduces T c . 

To be able to perform numerically demanding configurational averaging 
within the Dynamical Cluster Approximation (DCA), we use a model which 
excludes many realistic features. We assume a simple cubic lattice with a single 
s-state per site, which does not cover the triple degeneracy of impurity state of 



boron |l2j . The impurity potential is restricted to a single site having no long- 



range Coulomb tails, and eventual screening of thispotential is not assumed 



This model was already studied by Shirakawa et al [13[ within the BCS theory, 
using the Coherent Potential Approximation for disorder, and also studied by 
one of the present authors |14j in the DCA for uncorrelated boron distribution. 

The DCA allows us to include statistical weight of boron complexes given 
by their binary correlations. Nearest-neighbour positions of two boron atoms 
are favoured by a dimerization energy. To avoid large clusters, we assume that 
once the atom is in a dimer it cannot gain the dimerization energy from another 
neighbouring atom. For this model we evaluate the density of states on the 
Fermi level as a function of boron concentration and dimerization energy. The 
T c is then obtained from the McMillan formula. 

The paper is organized as follows. In section [2] we introduce the model 
Hamiltonian and describe the theory. In section [3] we calculate the density of 
states (DOS) and study the impact of boron correlations on No. In section [4] 
we compare our theoretical values of the critical temperature with experimental 
data. Section [5] contains conclusions. 

2. Theory 

The presence of boron impurities at concentrations reaching 5% of the crystal 
atoms modify many properties of the crystal ranging from the lattice constant, 
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over the phonon spectrum and electron-phonon coupling, to the dielectric func- 
tion. We neglect all these effects on superconductivity and express material 
parameters of McMillan formula as a function only of the density of states at 
the Fermi level No. We introduce a model Hamiltonian and a method to calcu- 
late the density of states. At the end of this section we specify the probability 
distribution of boron clusters. 

2.1. McMillan formula 

The McMillan formula 15] provides the critical temperature 
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Since the disorder-dependence of the Debye frequency in the boron-doped dia- 
mond is weak, we use the pure-diamond value Sxjd/^b = I860 K. 
The screened pseudopotential 



/i* = VN 



1 + VN In 



E c 



0.62 Hub 



(2) 



is given by the strength of the Coulomb interaction V. We have found neither 
measurement nor estimate of V in the literature, theqrefore we set this param- 
eter from experimental T c . A single value V = 3.6 x 10~ 23 cm 3 eV is used for all 
samples independently of boron concentration and correlation. 
We estimate the Coulomb cutoff energy as 



F X /V ( 9Ne 
Ec - 2 N ° {-dE 



E=0 



(3) 



For free particles this formula gives Ec — > h 2 kp/2m, usually used for the cutoff 
in metals 16]. In the impurity band the Fermi momentum k-p is not well- 
defined, but expression © in terms of the density of states is applicable. Since 
the Coulomb cutoff is only a subsidiary quantity, we use its uncorrelated value 
also for correlated samples. The concentration dependence must be maintained 
as one can see from the free-particle limit. 

The electron-phonon coupling A we take from ab initio calculations. The 
published results cover only a few concentrations, theqrefore we interpolate them 
in the spirit of the Morel- Anderson [17[ formula 



A = 



UN, 



(4) 



Here U represents a phonon-electron coupling strength and QNq describes a 
screening. 

Using formulae ([lJ-Sl) we obtain T c as a function of the density of states Nq 
and its derivative giving the Coulomb cutoff Eq . The Nq depends on the boron 
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concentration and dimerization energy. The Eq we take as a function of only 
the concentration. 



In [1J| it was shown that the Belitz theory with neglected vertex correc- 



tions provides a better agreement with data on 100 samples than the McMillan 
formula. For impurity concentrations up to 5% both approaches provide com- 
parable results. In this paper we thus use the McMillan formula which is the 
standard approximation. 



2.2. Density of states 

Now we specify a model from which we obtain the density of states of non- 
interacting electrons as a function of the boron concentration ne . It is given by 
a Hamiltonian of the valence band in diamond Hq and a random potential of 
boron impurities V, 

n = n + v = n Q + J2 r nSa]a u (5) 

i 

where r\i = 1 at impurity sites and zero elsewhere, and S is the potential ampli- 
tude. 

We employ the DCA for clusters of 2 x 2 x 2 atoms on a cubic lattice. 
We only briefly introduce the DCA here; all details can be found in a previous 
paper of one of the authors [3] or in the method's original paper 18]. The DCA 
provides a Green function averaged over all possible configurations of impurities 
on a cluster embedded in an effective medium. The key point is how to close 
the self-consistency constructing the effective medium from the averaged Green 
function. 

Naive constructions based on a tight-binding representation of the self-energy 
have turned out to violate analytic properties and thus causality. In the DCA 
the self-energy is approximated in the momentum representation. The Bril- 
louin zone is according to the size of the cluster divided into subzones - in our 
case 2x2x2 subzones. Within these subzones the self-energy is momentum- 
independent; that is, the selfenergy is represented by eight complex functions of 
frequency. 

All subzones contribute to the bare density of states 

K 

where K is a subzone index. In general, the subzone contribution p^r is obtained 
by integrating over the subzone with the electron dispersion of the valence band. 
For simplicity we approximate in each subzone K by a semielliptical function 

N c n w K y v u> K j 
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where N c is number of subzones and 

E K = k( E K* + E lt)> (°) 

reproduce the maximum Ej^ x and minimum E^ 11 energy in the subzonc K. 

The density of states © by definition maintains the width of the valence 
band. This approximation also yields a correct curvature at the edge; it correctly 
reproduces an effective mass of holes near the band edge. This feature is vital 
for a realistic description of the relatively shallow impurity state. 

The parameters wj^ and Ej^ are fitted to the tight-binding cosine band of 
width 22 eV. The impurity potential 6 = 8.91 eV reproduces the single-impurity 
bound state energy 0.37 eV above the valence band. 

Each boron atom releases a singe hole, thcqrcforc the density of holes equals 
n B . From this density we determine the Fermi energy E~p 

poo 

n B = 2 / p(E) dE. (10) 
Je f 

Here the factor of two accounts for spin. After we find E F we shift the energy 
reference point so that E-p = as it is usual in the theory of superconductivity. 

Mukuda et al Q argued that the concentration of holes can differ from the 
boron concentration due to boron atoms in neutral B-H complexes or interstitial 
positions. On the other hand Klein et al Q found that the effective number 
of carriers deduced from Hall-effect measurements was much larger than the 
number of boron atoms in samples. Since reliable hole concentrations are not 
accessible, we make the assumption that all samples are doped ideally. 

2.3. Boron correlations 

The exact form of boron correlations in the diamond remains a matter of 
discussion. We introduce their attractive correlation via an interaction energy 
£dim between boron atoms located at the nearest neighbour sites. 

Let I be a configuration of rix boron atoms on N sites of the cluster with 
qx boron pairs in nearest neighbour sites. Its probability 

p x = -^(l-ij^e^W (11) 
Z 

is given by a configuration statistical factor and by total energy scaled by the 
temperature of sample preparation T p . Typical T p for 100 and 111 samples 
is about HOOK. The probability is normalized to unity, Y)jPx — 1, which 
determines Z. The parameter x determines to the boron concentration 

n B = ^2 n iPi- ( 12 ) 
i 

For uncorrelated boron atoms, £dim — ^ 0, it equals the boron concentration, 
x — > n-Q. In our case N = 8 which allows for 2 8 configurations. By symmetry 
one can reduce a number of configurations which must be treated numerically. 
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Apparently, averaging over small clusters underestimates contribution of 
dimers. In the simplest case of 1 x 1 x 1 cluster the dimers are absent. In 
the assumed case of 2 x 2 x 2 all atoms are on the surface of the cluster with 
half of nearest neighbors in the cluster and half outside. The outer neighbors 
are not accounted for which reduces the effect of dimers crudely by half. It can 
be partly compensated by an artificial increase of the dimerization energy by 
£2 = -k B T p In 2. 

3. Density of states in the presence of boron dimers 

Energy levels of a boron dimer are rather distinct from a level of an isolated 
atom. While a hole binds to an isolated atom at 0.37 eV, for our model the 
dimer has the symmetric level at 1.6 eV and the anti-symmetric state forms only 
a resonant level in the valence band at —2.2 eV. For majority of samples the 
boron concentration is less than 5%, theqrefore isolated boron atoms prevail in 
uncorrelated configurations. When the energy gain of two boron atoms at neigh- 
bour positions exceeds the preparation temperature, Edim 3> K&T V ~ 94.8 meV, 
the fraction of dimers becomes large and this modifies the spectrum of energies 
in the sample. 

The effect of boron dimers on the overall density of states in the impurity 
band for 5% boron doping is shown in Fig. [TJ where we compare the DOS for 
uncorrelated case a) with two correlated cases. In the case b) we take the boron- 
boron binding energy £dim = — 160 me V resulting from ab initio calculations 
0, S fl(i| ]. We refer to this binding energy as a realistic correlation. In the 
case c) we use binding energy £dim = — 460 me V which yields the best fit of 
experimental data, see comparison in Fig. [5l We shal refer to this binding 
ebergy as a fitted correlation. 

It should be noted that ab initio studies compare energies of periodic crystals 
made of clusters B 2 C62 or B2C52 with selected configuration of boron atoms. 
Energies which one can interpret as the binding energy of a nearest neighbour 
dimer achieve values ranging from — 158 meV, over — 205 me V to — 288 me V 
depending on authors and eventual gradient corrections. The fitted correlation 
energy —460 meV is thus not far from these values. 

We focus on the impurity band because the Fermi energy, Ep = 0, lies there. 
The shape of the valence band is unimportant as it extends over energy interval 
from — 1.5eV to — 24.5eV. The thick line represents the density of states and 
thin lines are its decomposition into contributions according to number of boron 
atoms in the 2x2x2 cluster. All configurations of n boron atoms are added 
together so that they contribute with probability p n = *Y^x~ n Pi- 

Probabilities of n-boron clusters for the uncorrelated case are po — 0.66, 
Pl = 0.28 and P2 — 0.05. The zero-boron clusters contribute mainly to the 
valence band. In the impurity band they give a small contribution due to states 
tunneling from the surrounding medium. States from the single-boron clusters 
dominate the impurity band. 

The two-boron clusters include 12 nearest-neighbour configurations which 
are the dimers, 12 from second-nearest neighbours, and four from third-nearest 
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Figure 1: Density of states at doping ?ib = 0.05. Positions of boron atoms at neighbour sites 
are a) uncorrelated b) correlated by energy £dim = — 160 me V found from realistic ab initio 
studies and c) correlated by e<jim = —460 meV fitted to give best agreement with experimental 
data. The growth temperature is T p = 1100 K. Thick line represents the DOS and thin lines 
its decomposition into contributions according to number of boron atoms in the 2x2x2 
cluster. ' 



neighbours. For uncorrelated system in Fig. QJi) the dimers increase the DOS 
at the upper edge of the impurity band around E ~ 1.55 eV, and in the va- 
lence band at E ~ —2.15 eV. In the correlated cases in Fig. [Ud) and c) their 
contribution is even more pronounced. The position of the two-boron levels are 
shifted due to an increase of the dimer fraction in two-boron clusters and due 
to feedback effect of surrounding medium. 

The small split band near E = 2.4 eV in uncorrelated case in Fig. QJi) is 
due to three-boron complexes. For realistic correlations it shifts towards higher 
energies. For fitted correlations it merges with the main impurity band. States 
of clusters with more boron atoms are also present, but their contributions are 
very small for the given concentration. 

The realistic correlated case b) differs from the uncorrelated one a) mainly 
in the probability of zero-, single-, and double-boron cluster contributions. For 
given binding energy £di m = — 160 me V they are po — 0.70, p\ — 0.21 and 
P2 = 0.08. In accordance with decreased p\ and increased p%, the contributions 
of double and higher boron states are enhanced and the peak associated with 
single boron states is reduced. The lower probability of the single-boron clusters 
leads to the decrease of the DOS at the Fermi energy. 

With fitted correlation c) the contributions from larger clusters of boron 
atoms become more apparent. For £<jim = — 460 me V we find po — 0.77, p\ 
0.07 and P2 — 0.15, theqrefore the number of boron atoms in dimers is four 
times higher than the number of isolated atoms. As one can see in Fig. [TJ;) the 
major part of the impurity band is due to dimer states. 

Note that even when dimers dominate, the Fermi energy remains in the 
single-atom part of the impurity band. This follows from the fact that each 
dimer binds two holes of spin up and down in the symmetric state. Since the 
number of boron atoms in dimers equals the number of bounded holes, two 
boron donors become passivated when they form a dimer. The metallic band is 
thus formed only by unpaired boron atoms. 

The dimers influence the density of states in two ways. First, they effectively 
reduce the doping level. Second, they contribute to the medium in which the 
band is formed. To demonstrate the second effect in Fig. [5] we have included a 
DOS for uncorrelated system with ub = 0.0091 giving the probability of single- 
atom clusters p\ — 0.07 being equal to p\ of the fitted correlation. One can 
see that the approximation of a correlated system by a non-correlated one with 
dimers subtracted from the boron concentration leads to correct shape of the 
metallic band except for the contribution of clusters with dimers to the density 
of states at the Fermi level. For parameters in Fig. [5] this contribution is a 
significant ~ 20%, assuming exponential dependence of T c on N . 

The dependence of N on the strength of boron correlations is shown in 
Fig. [3] One can see that positive correlations in boron distribution lower the 
density of states at the Fermi level N for any doping Fig. [3] also presents 
the Coulomb cutoff given by formula © . This value is very close to the distance 
of the Fermi level from the right band edge. An alternative definition of the 
Coulomb cutoff from the bad edge, however, is imposed due to split off bands 
of trimers and larger clusters of low probability. 
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Figure 2: Approximation by subtracted dimers. The density of states at doping no = 0.05 
and fitted correlation £dim = —460 meV (gray re-plot from Fig. [TJ:) is compared with uncor- 
rected system of concentration = 0.0091 (in black lines). Both systems have the same 
probability of single-boron clusters giving similar 1DOS (thin line). They have also similar 
0DOS (dotted lines). The single-boron states, however, easily tunnel into two-boron clusters 
giving appreciably larger contribution of 2DOS (dashed lines) at the metallic band. 




Figure 3: The density of states at the Fermi energy No as a function of boron concentration 
«B f° r no (f un line), realistic (dashed-dot line) and fitted (dashed line) correlation. The left 
vertical axis corresponds to the Coulomb cutoff Eq (dotted line) evaluated for uncorrelated 
boron atoms. 
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Figure 4: Coupling constant A as a function of the boron concentration. Symbols represent 
results of ab initio calculations of A for uncorrelated boron-doped diamond. Solid line re- 
produces A dependence on no using formula ||4} with uncorrelated No(n-g) (see Fig. [3J and 
parameters Q = 35.7 X 10 — 23 cm 3 eV, U = 23.3 X 10 — 23 cm 3 eV. Ab initio results were obtained 
within the Virtual Crystal (+ [H], X HJ, * [H) and Super Cell (□ 0, A ||, O 0) 
method. 

4. Dependence of critical temperature on boron correlations 

Now we are ready to discuss the impact of correlations on T c . We ask the 
question whether different T c of the 100 and 111 samples can be explained by 
a different content of dimers. Mukuda et al [8[ have observed a strong NMR 
line of single-boron state in the 111 samples while this line was rather weak 
in the 100 samples. Based on this experimental fact we treat 111 samples as 
uncorrelated while in the 100 and HPHT samples we include dimers. 

First we specify material parameters. The Q and U are chosen to re- 
produce ab initio calculations of coupling parameter A in uncorrelated boron- 
doped diamond. Formula ((3]) with uncorrelated N (n^) and parameters Q = 
35.7 x 10 _23 cm 3 eV and U = 23.3 x 10 _23 cm 3 eV is compared with ab initio values 
in Fig. 2] Values V and £di m are not accessible in print, we thus fit them to exper- 
imental values of T c . The best agreement is achieved for V = 3.6 x 10~ 23 cm 3 eV 
and £dim = —460 meV. 

In Fig. [5] we compare T c from McMillan formula ([1} for uncorrelated and 
correlated systems with experimental data for 111, 100 and HPHT samples. 
As one can see the realistic correlation energy £dim = — 160 me V found in ab 
initio studies is too low to explain differences in the critical temperature. With 
the fitted correlation energy £dim = —460 meV, the theory reproduces trends 
of experimental values for small concentrations. For iib > 0.02 one finds only 
qualitative agreement. 

Compared to the experimental data, the present study underestimates the 
difference between 111 and 100 samples. It can be in part due to our simple 
model, in part due to approximations employed to evaluate the density of states. 
First of all, inside the 2x2x2 cluster each atom has only three nearest neigh- 
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Figure 5: Critical temperature T c as a function of boron doping for no (solid line), realistic 
(dashed-dot line), and fitted correlation (dashed line). Open symbols represent 111 MPCVD 
data (□ H, 0, O 0. A [H, v 01), full symbols 100 MPCVD data (. [J], A @, ■ pj, ♦ 
[13, ▼ [H) and crosses HPHT samples (+ [l[, X [H, * [U). 



bours while in the crystal it has six. Accodingly, each boron can form only three 
dimers inside the cluster while three dimers cross the cluster border. The cor- 
relations corresponding to the former are covered in our treatment while those 
corresponding to the latter are not. 

There are also additional mechanisms that appear already in our model but 
have not been included for simplicity of discussion. For example, two boron 
atoms at second-neighbour positions form a symmetric state of energy 0.91 eV 
which is close to the single-atom bound state and at concentrations above 1% 
contributes to the mettalic band. The anti-symmetric state of such a boron 
pair has energy at — 0.8 eV which forms a resonant level in the valence band 
and does not contribute to the impurity band. One can crudely say one half of 
boron atoms in the second-neighbour dimers are passivated. 



5. Conclusion 

In this paper we have addressed the impact of correlations of boron impu- 
rities on the density of states and thus on the critical temperature T c of the 
superconducting transition. To this end we have calculated the impurity-band 
density of states for several strengths of attractive correlations between boron 
atoms. Calculations have been done within a one-band cubic lattice with on-site 
impurity model employing the DCA method on a cluster of 2 x 2 x 2 atoms. 

Wc have found that correlations enhance states associated with multi-boron 
complexes reducing the density of states at the Fermi energy, see Fig. Q] and 
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Fig- El Following NMR data we have approximated 111 oriented samples as an 
uncorrelated system and 100 samples as a system with correlations. We have 
found that an ab initio binding energy £di m = — 160 me V of boron dimer is 
not sufficient to explain differences in T c . From the best fit we have found the 
binding energy £di m = —460 meV which explains experimental data at least for 
boron concentrations below 2%. 

The studied model has several shortcomings. It has overly-simplified elec- 
tronic band structure and the binding potential is single-site. Also our treatment 
of the disorder covers only contributions of boron complexes to the density of 
states but neglects other disorder effects like the weak localization, which has 



been discussed as a possible mechanism enhancing the superconductivity. [31 [ 
Nevertheless, our results strongly indicate that binary correlations are respon- 
sible for different critical temperatures of the 111 and 100 samples. 
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